Microscopic Structural Relaxation in a Sheared Supercooled Colloidal Liquid 
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The rheology of dense amorphous materials under large shear strain is not fully understood, partly 
due to the difficulty of directly viewing the microscopic details of such materials. We use a colloidal 
suspension to simulate amorphous materials, and study the shear-induced structural relaxation with 
fast confocal microscopy. We quantify the plastic rearrangements of the particles in several ways. 
Each of these measures of plasticity reveals spatially heterogeneous dynamics, with localized regions 
where many particles are strongly rearranging by these measures. We examine the shapes of these 
regions and find them to be essentially isotropic, with no alignment in any particular direction. 
Furthermore, individual particles are equally likely to move in any direction, other than the overall 
bias imposed by the strain. 

PACS numbers: 82.70.Dd, 61.43.Fs, 83.60.Rs 



I. INTRODUCTION 

Many common materials have an amorphous struc- 
ture, such as shaving cream, ketchup, toothpaste, gels, 
and window glass [11 0, S I3|- -"- n some situations these 
are viscous liquids, for example when window glass is 
heated above the glass transition temperature, or a shav- 
ing cream foam that has been diluted by water to become 
a liquid with bubbles in it. In other situations these are 
viscoelastic or elastic solids, such as gels and solid win- 
dow glass [j| • For solid-like behavior, when a small stress 
is applied, the materials maintain their own shapes; at 
larger stresses above the yield stress, they will start to 
flow @, 0, Hi • Understanding how these materials yield 
and flow is important for the processing of these mate- 
rials, and understanding their strength in the solid state 

A particularly interesting system to study is a colloidal 
suspension. These consist of micron or submicron sized 
solid particles in a liquid. At high particle concentra- 
tion, macroscopically these are pastes and thus of practi- 
cal relevance [6J. Additionally, for particles with sim- 
ple hard-sphere like interactions, colloidal suspensions 
also serve as useful model systems of liquids, crystals, 
and glasses [H, HH, [TJ, Il5l |. Such colloidal model sys- 
tems have the advanta ge t hat they can be directly ob- 
served with microscopy particular inter- 
est in this paper is using colloidal suspensions to model 
supercooled and glassy materials. The control param- 
eter for hard sphere systems is the concentration, ex- 
pressed as the volume fraction <fi, and the system acts 
like a glass for <j) > 4>g ~ 0.58. The transition is the 
point where particles no longer diffuse through the sam- 
ple; for cf> < <j) g spheres do diffuse at long times, al- 
though the asymptotic diffusion coefficient decreases 
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sharply as the concentration increases [19|, [2Cj, |21| . The 
transition at </> s occurs even though the spheres are not 
completely packed together; in fact, the density must 
be increased to iArpp » 0.64 for "random-close-packed" 
spheres [22|, [23|, \24, |25|. [26| before the spheres are mo- 
tionless. Prior work has shown remarkable similarities 
between colloidal susp ensions and conventional molecu- 
lar glasses 0, m, he m, m, h m, m im . 

One important unsolved problem related to amorphous 
materials is to understand the origin of their unique rhe- 
ological behavior under shear flow. Early in the 70's, 
theory predicted the existence of "flow defects" beyond 
yieldin g [9|1 , later termed shear transformation zones 
(STZ) 34]. These microscopic motions result in plas- 
tic deformation of the sheared samples [H, [36[. Simu- 
lations later found STZs by examining the microscopic 
local particle motions (TT1 . 137L l38j ] . Recently fast confo- 
cal microscopy has been used to examine the shear of 
colloidal suspensions [39| . and STZ's have been directly 
observed [lj3]. This provided direct evidence to support 
theoretical work on STZ's @, O H3, 53 ■ 

However, questions still remain. First, most of the 
prior work has focused on the densest possible sam ples , 
at concentrations which are glassy (<f) > 4>g) fl5l. |39|. 
Given that the macroscopic viscosity of colloidal suspen- 
sions change dramatically near and above <f>c [421 ] . it is 
of interest to study slightly less dense suspensions un- 
der shear, for which rearrangements might be easier [43| . 
In this paper, we present such results. Second, prior in- 
vestigations of sheared amorphous materials have used 
a va riety of different ways to quantify plastic deforma- 
tion [TBI. [39L l44l [4~5l] . In this paper, we will compare and 
contrast plastic deformations defined in several different 
ways. While they do capture different aspects of plastic 
deformation, we find that some results are universal. In 
particular, in a sheared suspension, there are three dis- 
tinct directions: the strain velocity, the velocity gradient, 
and the direction mutually perpendicular to the first two 
(the "vorticity" direction). We find that plastic deforma- 
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tions are isotropic with respect to these three directions, 
apart from the trivial anisotropy due to the velocity gra- 
dient. The deformations are both isotropic in the sense 
of individual particle motions, and in the sense of the 
shape of regions of rearranging particles. 



II. EXPERIMENTAL METHODS 

The experimental setup of our shear cell is shown in 
Fig. [T]and is similar to that described in Refs. [39|, |46| |. 
The glass plates are 15 mm in diameter, and to the top 
plate is glued a small piece of glass with dimensions 
5 mmxl mm, with the long dimension oriented in the 
direction of motion x. The purpose of this piece of glass 
is to decrease the effective gap size. Between the plates 
are three ball bearings, used to control the gap size; for 
all of our data, we maintain a gap size of H = 130 /im. 
Over the 5 mm length of the small pieces of glass, the 
gap varies by no more than 15 /xm; over the narrower di- 
mension of the small pieces of glass, the gap varies by no 
more than 10 /urn. Thus, overall the sample is between 
two plates which are parallel to within 1% and in the 
direction of shear, they are parallel to within 0.3%. 

A droplet of the sample (volume ^200 fA) is placed 
between the two pieces of glass. The top plate is free 
to move in the x direction, and the bottom plate is mo- 
tionless. The shear rate is controlled by a piezo-electric 
actuator (Piezomechanik GmbH Co.) driven by a trian- 
gular wave signal with a period ranging from T = 150 
to 450 s and an amplitude of A = 175 /im. Thus we 
achieve strains of 70 = A/H = 1.4. Prior to taking data, 
we allow the shear cell to go through at least one com- 
plete period, but usually not more than three complete 
periods. 
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FIG. 1: (Color online) Sketch of the shear cell. A fluorescent 
sample (gray, green online) is put between two parallel glass 
plates (dark gray, blue online) with gap H set by three ball 
bearings, two of which are shown. The top plate is movable, 
controlled by a piezo motor driven by a triangular wave. The 
bottom plate is fixed, and the confocal microscope takes im- 
ages from underneath. Note the definition of the coordinate 
system, where x is in the velocity direction and z is in the 
velocity gradient direction. 

Our samples (Fig. [J) are poly(methyl methacrylate) 
colloids sterically stabilized with poly-12-hydroxystearic 
acid [47| | . These particles are suspended in a mixture 
of 85% cyclohexybromide and 15% decalin by weight. 
This mixture matches both the density and the index of 



refraction of the particles. To visualize the particles, they 
are dyed with rhodamine 6G [48| . The particles have a 
radius a = 1.05 ± 0.04 /im, with the error bar indicating 
the uncertainty in the mean diameter; additionally the 
particles have a polydispersity of no more than 5% (as 
these particles can crystallize fairly easily) [T^.l49l. [5fl[5l1 |. 

In this work we study several samples with volume 
fractions between 0.51 and 0.57. Thus, our samples 
are quite dense liquids, comparable to prior work with 
"supercooled" colloidal liquids [12, HH • The differences 
in volume fraction between samples are certain to within 
±0.01, and the absolute volume fraction has a system- 
atic uncertainty of ±0.06 due to the uncertainty of the 
particle radius a. However, none of our samples appear 
glassy, and thus we are confident our maximum volume 
fraction is less than ~ 0.6. While the particles in decalin 
behave as hard-spheres, in our solvent mixture they carry 
a slight charge. This does not seem to affect the phase 
behavior dramatically at high volume fractions such as 
what we consider in this work; see for example |53| . 

To characterize the relative importance of Brownian 
motion and the imposed strain field, we can compute the 
modified Peclet number, Pe* — r ya 2 /2D OQ , where is 
the long time diffusion coefficient of the quiescent sam- 
ple. We measure from mean square displacement 
data taken from the quiescent sample with the same vol- 
ume fraction. The large At data for the mean square 
displacement can be fit using (Ax 2 ) = 2D oc At. Roughly, 
a 2 /2Z3 oc reflects the average duration a particle is caged 
by its neighbors in the dense suspension. 

7 is the imposed strain rate 2A/(HT), and a is the 
particle size. The extra factor of 2 in 7 is because we use 
a triangle wave, and thus the half period sets the strain 
rate. For our samples, we find Doo w 5 x 10~ 4 /im 2 /s, 
and we have 7 ranging from 0.0060 to 0.0180 s _1 ; thus 
Pe* « 7 — 20. Given that Pe* > 1, the implication is 
that the motions we will observe are primarily caused by 
the strain, rather than due to Brownian motion. We use 
the modified Peclet number based on rather than the 
bare Peclet number based on the dilute-limit diffusivity 
Do, as we will focus our attention on the dynamics at long 
time scales, which we will show are indeed shear-induced. 

Shear-induced crystallization has been found in previ- 
ous work [H, As we wish to focus on the case of 
sheared amorphous materials, we check our data to look 
for crystalline reg ions using standard methods which de- 
tect ordering [52l l56l. l57l, 58]. Using these methods, 
we find that particles in apparently crystalline regions 
comprise less than 3% of the particles in each of our ex- 
periments, and are not clustered, suggesting that the ap- 
parently crystalline regions are tiny. This confirms that 
our samples maintain amorphous structure over the time 
scale of our experiments, although perhaps if we contin- 
ued the shearing over many more cycles we would find 
shear-induced crystallization. 

We use a confocal microscope to image our sample (the 
"VT-Eye", Visitech), using a 100 x oil lens (numerical 
aperture = 1.40) [H |H, H|. A 3D image with a vol- 



3 




FIG. 2: A 20x20x20 /im 3 image of a supercooled colloidal 
liquid taken 15 /im away from the fixed bottom plate by our 
confocal microscope in less than 1 s. The scale bar represents 
10 /mi. 



ume 50x50x20 /im 3 is acquired in less than 2 s; these 
images contain about 6000 particles. The 3D image is 
256 x 256 x 100 pixels, so approximately 0.2 /im per pixel 
in each direction. Figure [2] shows a representative image 
from a somewhat smaller volume. The 2 s acquisition 
time is several orders of magnitude faster than the diffu- 
sion for particles in our high volume fraction sample. To 
avoid any boundary effects (fjoj , we scan a volume at least 
20 /im away from the bottom plate. Particle positions are 
determined with an accuracy of 0.05 /im in x and y, and 
0.1 jum in z. This is done by first spatially filtering the 
3D image with a bandpass filter designed to remove noise 
at high and low spatial frequencies, and then looking for 
local maxima in the image intensity 16 i f . Our tracking 
algorithm is similar to prior work |48t l6l[ |62| , where we 
first identify particles within each 3D image, next remove 
the overall average shear-induced motion from all of the 
particles, then track the particles in the "co-shearing" 
reference frame using conventional techniques [6lj . and 
finally add back in the shear-induced motion that was 
previously removed. This is similar to the "Iterated CG 
tracking" method described in Ref. [62]. The key idea 
of this tracking is that particles should not move more 
than an inter-particle spacing between each image; this 
condition is satisfied in the co-shearing reference frame. 

Due to the strain, particles that start near one face of 
the imaging volume are carried outside the field of view, 
while on the opposite face new particles are brought in- 
side. Thus, for larger strains, the total number of par- 
ticles viewed for the entire duration diminishes. For the 
data discussed in this work, we consider both instanta- 
neous quantities and quantities averaged over the entire 



half-cycle of strain. For the former, we view ~ 5500 par- 
ticles, while for the latter, we typically can follow ~ 3000 
particles, which limits our statistics slightly. 



III. RESULTS 
A. Locally observed strain 

Our goal is to understand if the local shear-induced 
motion is isotropic in character. However, first we seek 
to understand and quantify the more global response of 
our sheared samples. 

When shearing disordered materials or complex flu- 
ids, one often finds shear localization or shear banding, 
due to the nonlinear yielding and relaxation in local re- 
gions [fj, lZ|, l63j, l64|. To check for this in our data, we 
start by taking 3D images under the applied shear rate 
7macro = 0.016 s _1 over a very large range in z, from 
to 70 /im away from the bottom plate, almost half of 
the gap between two shearing plates. To allow us to vi- 
sualize more clearly over such a large depth, we dilute 
the dye concentration by mixing dyed and undyed col- 
loids at a number ratio of around 1:80 and keeping the 
desired volume fraction <f> « 0.50. Our sample does in- 
deed form a shear band, as shown in Fig. [31 which shows 
the particle velocity v x in the direction of the shear as 
a function of the depth z. The velocity changes rapidly 
with z in the range < z < 20 /im, and then more 
slowly for z > 20 /im; thus much of the shear occurs 
adjacent to the stationary plate at z — /im, similar 
to prior work which found shear adjacent to one of the 
walls d, 0, [H, [H, [6f| . Furthermore, the velocity profile 
is relatively stable during the course of the half-period, 
as seen by the agreement between the velocity profiles 
taken at different times during this half period (different 
lines in Fig. [3] Thus, the shear band develops quickly in- 
side the supercooled colloidal liquid, and remains fairly 
steady under the constant applied strain rate. The loca- 
tion and size of the shear band varies from experiment 
to experiment. 

Given the existence of a shear band, the applied strain 
is not always the local strain. In this work we wish to 
focus on the motion induced by a local strain, rather 
than the global formation of shear bands. Thus, for all 
data sets presented below, we always calculate the local 
instantaneous strain rate j me so(t) = V;c ^ z+ ^ z ^>~ Va: ^ z ^ , 
Here Az = 20 /zm is the height of the imaged volume. 
Related to 7 m0 so> we can calculate the total local applied 
strain by integrating j m eso{t): 

7mcso(£) = / 7meso(*') d *'- W 

Jo 

Furthermore, we verify that for each data set considered 
below, v x (z) varies linearly with z within the experimen- 
tal uncertainty, and thus 7 meS o is well-defined, even if 
globally it varies (Fig. [3]) . 
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FIG. 3: (Color online) Shear profile v x (z) within a 
50x50x70 /im 3 volume, measured at different times during 
one half period T/2 — 84 s. The total macroscopic strain 
applied in 84 s is 1.35. The sample has a volume fraction (j> m 
0.50, and the global applied shear rate is 7macro = 0.016 s _1 
controlled by the piezo motor. If there was not a shear band, 
one would expect to observe a linear velocity profile as in- 
dicated by the lower thick diagonal line. The four curves 
represent different times during the half period as indicated. 
Note that these data are obtained using a coarse-grained im- 
age velocimetry method [fiH], rather than particle tracking. 
Here we take rapid images over the full volume with a spac- 
ing of z = 1 fj,m in the vertical direction. For each value of 
z, we cross-correlate subsequent images to obtain a mean in- 
stantaneous velocity v x with a resolution of 0.05 /im/s, set by 
the pixel size and the time between images. See Ref. [621 ] for 
further details. 



paper, the choice t = T/2 will be assumed, except where 
noted when we wish to characterize the mesoscopic strain 
for time scales shorter than T/2. 




*(S) 

FIG. 4: A typical example of a measured local instanta- 
neous shear rate over one period of shearing in a gap of 
20—40 fim away from the fixed bottom plate. The sam- 
ple has volume fraction <f> = 0.51, the applied strain rate is 
Tmacro = 0.013 s _1 , and the period is T = 200 s. 

At this point, we have defined the key control param- 
eters, which are measured from each experimental data 
set: the strain rate 7 meS o and the strain amplitude 7 mcS o- 
We next consider how these variables relate to the mag- 
nitude of the shear-induced particle motion. 



While the shape of the shear band is essentially con- 
stant (Fig. [3|), in many cases the local strain rate varies 
slightly with time. As our forcing is a triangle wave, over 
any given half cycle the global applied strain rate is a 
constant. We can measure the local strain rate for each 
data set; a typical example is shown in Fig. 2J This fig- 
ure shows the local instantaneous measured strain rate 
within the region 20 /im < z < 40 /zm, over one full 
period. After the shearing starts, 7 mcS o quickly rises up 
to 0.015 s _1 , implying that the shear band has formed. 
Subsequently, the local strain rate continues to increase 
up to 0.030 s" 1 over the rest of the half period. The 
small fluctuations of 7meso are due to the microscopic re- 
arrangements of particles, which can be somewhat inter- 
mittent. Given that the local instantaneous strain rate is 
not constant (despite the constant applied strain rate), 
we will characterize our data sets by the time averaged 
local strain rate defined as 



typically using t = T/2, the half period of the strain. 
That is, we consider 7 mC so as one key parameter charac- 
terizing each data set, although we will show that we see 
little dependence on this parameter. In the rest of the 



B. Individual particle motions 

Because of our large local strains (measured to be 
7mcso > 0.3 for all cases), we observe significant particle 
motion, as shown Fig.(5Ja,b). In the laboratory reference 
frame, the microscopic velocity gradient is obvious either 
in the raw trajectories [Fig. 02a)] or in the large displace- 
ments [Fig. EJb)] measured between the beginning and 
end of the half period. However, in a sense, much of this 
motion is "trivial"; we wish to observe what nontrivial 
local rearrangements are caused by the strain. To do this, 
we consider the non-affine motion by removing averaged 
particle displacements at the same depth z from the real 
trajectories of particles [38|, [39|, |4J, |66|, [67( , 

Xi{t)=Xi{t)- [ j meso (t')zi(t')dt' (3) 
Jo 

where the removed integral represents the shearing his- 
tory of the particle i. To be clear, the shearing history is 
based on the average motion within the entire imaged re- 
gion, and the remaining motion of particle i is caused by 
interactions with neighboring particles due to the shear. 
This motion is shown in Fig. EJc) , showing the xz plane 
rather than the xz plane; the particles move shorter dis- 
tances. Their overall displacements in this "de-sheared" 
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reference frame are shown in Fig.^d). A few trajectories 
are up to 2 /im long, comparable to the particle diame- 
ter. These non-affine displacements shown in Fig. [5](d) 
are much smaller than the raw displacements of Fig.[5fb), 
but much larger than thermally activated Brownian mo- 
tion, which takes more than 1000 s to diffuse over a 1 /im 
distance in our dense samples (comparable to the parti- 
cle radius a). These non-affine motions x reflect shear- 
induced plastic changes inside the structure. 



of shear. Figure EJa) shows mean squared displacement 



(MSD) of the non-affine motion Af = Aa 



Ay 2 




FIG. 5: Trajectories of colloids in an xz slice (5 /im thick 
in the y direction), see Fig. [1] for the coordinate axes. The 
sample has (j> — 0.51, and the data shown correspond to a 
locally measured accumulated strain 7 meS o = 0.43 over 45 s 
of data, so the effective strain rate is 7 moS o = 7 mcS o/At = 
0.0096 s" 1 . (a) Trajectories in a reference frame co-moving 
with the average velocity v x (z — 20 /im). (b) Displacements 
corresponding to data from panel (a), where the start point 
of each particle is marked with a circle and the end point is 
marked with an arrowhead, (c) The same data, but with the 
afSne motion removed; this is the xz plane, (d) Displacements 
corresponding to panel (c). 

To quantify the amount of this non-affine motion x, one 
could calculate the mean squared displacements (MSD) 
often defined as 



(Ai 2 )(At) = ((xi(t + At) 



(4) 



where the angle brackets indicate an average over time t, 
as well as particles i. Thus this identity assumes that 
environmental conditions remain the same for all the 
time, since it does not depend on t. However, as shown 
by Fig. U the shear rate 7iocai(i) depends on the time. 
Therefore, we use an alternate formulation 



(Ax')(t) = ((x i (t)-x i (t = 0)Y) i 



(5) 



where the angle brackets only indicate an average over 
particles, and t is the time since the start of a half period 



Az as a function of t for five different experiments with 
different strain rates 7 mcS o, from 0.02 to 0.006 s -1 . In 
each case the curves nearly reach a slope of 1 on the log- 
log plot, indicating that shear quickly facilitates particles' 
rearrangements. The magnitude of the motion is Af 2 « 
1 /xm 2 , indicating that the original structure is mostly 
lost [681 ]. 

Figure [S^a) also shows that Af 2 is larger for faster 
strain rates at the same t. We find that this motion 
is determined by the accumulated strain, as shown in 
Fig. (H^b) by replotting the MSD as a function of 7 mcS o 
(Eqn. [1]). In this graph, the curves are grouped closer 
together and there is no obvious dependence on 7 mcS o- 

It suggests that the accumulated strain is an important 
parameter in the structural relaxation, which was also 
found in previous work on shear transformation zones 
[67l |69| , and is similar behavior to that seen for athermal 
sheared systems [7(j. Additionally, Fig. [Sfb) shows that 
the slopes of the curves are close to 1 when 7 mcS o > 0.1, 
confirming that the accumulated strain in our experi- 
ments is large enough to rearrange the original struc- 
ture in a supercooled colloidal liquid. We stress that the 
rough agreement between the curves seen in Fig. [5{b) is 
based on the locally measured applied strain, and not the 
macroscopically applied strain. 

An earlier study of steady shear applied to colloidal 
glasses by Besseling et al. [39| found that the diffusion 
time scale r scaled as 7 -0 ' 8 , and simulations also found 
power law scaling [53, HU. The collapse of our MSD 
curves [Fig. Ob)] seems to imply r ~ 7 -10 . It is pos- 
sible that the disagreement between these results is too 
slight to be clear over our limited range in 7 (less than 
one decade). Also, we study supercooled fluids whereas 
Ref. [1^] examines colloidal glasses. Furthermore, our 
maximum local strain is 7 mC so = 1-6, while Ref. (39j con- 
siders steady strain up to a total accumulated strain of 
10. Another recent study of sheared colloidal glasses [Til] 
implies a result similar to ours, r ~ 7 "_j but did not 
discuss the apparent difference with Ref. [391 ] . 

To better understand the mean square displacement 
curves (Fig. [6]), we wish to examine the data being aver- 
aged to make these curves. We do so by plotting the dis- 
tributions of displacements Af in Fig. [7J To better com- 
pare the shapes of these distributions, we normalize these 
displacements by the strain and thus plot P(As) where 
As = Af/ (7moso)' this normalization is motivated by the 
observation that at large 7 meS o, we have (Af 2 ) ~ 7mcso 
(Fig. [5]) [67|. Furthermore, we normalize P(As) so that 
the integral over all vectors As is equal to 1, similar to 
Ref. [67| . Figure [7J shows that the distributions corre- 
sponding to small strains are much broader than those 
corresponding to large strains. For the smallest strain, 
the distribution has a large exponential tail over 3 orders 
of magnitude. For larger strains (7 mcso > 0.1), the curves 
are no longer exponential and the tails are shorter, indi- 
cating fewer extreme events. These curves appear more 
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FIG. 6: (Color online) (a) The mean square displacement of 
the non-affine motion Af as a function of the time t since 
the start of shear, (b) The same data plotted as a function 
of the accumulated strain 7 meS o = / A /mcso(t')dt' . The five 
curves represents five data sets with the same volume fraction 
4> = 0.51 but three different shear rates 7 moS o as indicated in 
the figure. The total accumulated strain is the final point 
reached by each curve in panel (b). 



like Gaussian distributions. At the larger strain values 
(7mcso > 0.1), the distributions collapse; this is the same 
strain regime for which the mean square displacement 
becomes linear with 7 meS o [Fig- E^b)]. As the Ref. [67} 
suggests, the exponential tail for small strains is simi- 
lar to what has been seen for individual plastic events 
[6^ . |72| , while the distributions for larger strains are con- 
sistent with successive temporally uncorrelated plastic 
events drawn from the exponential distribution. How- 
ever, it is possible that these events are spatially corre- 
lated, which will be seen below. 

The mean square displacement data we have shown 
(Fig. [5]) treats all three directions equivalently, with the 
exception that the x displacements have had their non- 
affine motions removed. However, the three directions 
are not equivalent physically: the x direction corresponds 
to the shear-induced velocity, the y direction is the vor- 
ticity direction, and z is the velocity gradient direction. 
To look for differences in motion between these three di- 
rections, we plot the probability distribution of the dis- 
placements [Ax, Ay, Az], in Fig. [5) The three distribu- 
tions agree with each other, and in fact are symmetric 
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7: (Color online) Probability distribution function of 
Af/7me SO , the nonafBne displacement Af scaled by the 
local strain 7mcso- The data shown are from one experimental 
run, using portions of the data corresponding to strain incre- 
ments 7mcso = 0.01, 0.02, 0.10, and 0.43. These correspond 
to time intervals At = 1.5, 4.5, 16.5, and 45 s. P(As) is nor- 
malized so that the integral over all vectors As is equal to 1, 
similar to Ref. [67l ]. For small strains, the curves have large ex- 
ponential tails, and for larger strain, the tails become smaller 
and appear more Gaussian. The dashed line is Gaussian fit 
for 7meso = 0.43. For this sample, the data and parameters 
are the same as Fig. [5] 



around the origin. This suggests that the shear-induced 
motions are isotropic. Furthermore, they are well-fit by a 
Gaussian, suggesting that the shear-induced motion liq- 
uefies the sample (at least at the large 7 meS o considered 
for Fig. [5]). This seems natural in the context of jamming, 
where adding more strain moves the sample farther from 
the jammed state @. Of course, in our raw data the Ax 
data show a significant bias in the direction of the shear- 
induced velocity; but it is striking that the non-affine 
displacements Ax show no difference from the displace- 
ments in y and z, as also found in sheared colloidal glass 

Thus far, we have established that shear-induced par- 
ticle displacements are closely tied to the total applied 
strain 7 meS o [Fig.H^b)]. We then introduced the nonaffinc 
motion Af which we find to be isotropic, on the particle 
scale: individual particles are equally likely to have shear- 
induced displacements in any direction (Fig. [5]). While 
the distributions of displacements are isotropic, this does 
not imply that displacements are uncorrelated spatially. 
To check for this, we calculate two displacement correla- 
tion functions as defined in Ref. 173, U3, LZa : 



S P (R, At) 
S Sr (R, At) 



(Afj-Afj) 
(Af 2 ) ' 

asm ' 



(6) 
(7) 



where the angle brackets indicate an average over all par- 
ticles; see Refs. [73|, [75| for more details. The mobility is 
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Ax, A j, Az (|xm) 

FIG. 8: (Color online) Probability distribution functions for 
non-affine motions in each direction, Ax (filled cycle), Ay 
(hollow circle) and Az (hollow triangle). The red curve is a 
Gaussian fit to the Ax data. For this sample, the parameters 
are the same as Fig. \E\ (tj> = 0.51, 7 mcso = 0.43, 7 mcso = 
0.0096 s"\ At = 45 s). 




FIG. 9: (Color online) Spatial correlation functions Sf and 
Ssr characterizing particle motion. The data are the same 
as Fig. \5\ (4> — 0.51, 7 mcS o = 0.43 over 45 s of data, and 
Tmeso = 7mcso/At = 0.0096 s _1 ). The solid and dashed lines 
correspond to Sf and Ss r averaged over all the particles in 
the sample. The dotted lines (red) correspond to exponential 
fits with length scales £p = 6.4 /im and £$ r = 2.5 /im. 



that neighboring particles are correlated in their mo- 
tion. In particular, the vector-based correlation S? has 
a large magnitude at small R, showing neighboring par- 
ticles have strongly correlated directions of motion, even 
given that we are only considering the nonaffine displace- 
ments. The two correlation functions decay somewhat 
exponentially, as indicated by the straight line fits shown 
in Fig. [5J with decay constants £p = 6.4 /jm= 6.1a and 
£sr = 2.5 fim= 2.4a (in terms of the particle radius a). 
The larger the slope, the more localized the correlation 
is. is similar to that found previously for supercooled 
colloidal liquids, and £$ r is slightly shorter than the prior 
results 73]. Overall, these results confirm that the shear- 
induced particle motion is spatially heterogeneous, quite 
similar to what has been seen in unsheared dense liquids 
[H M, M, M, llZl H Izl and granular materials [HJ HI • 
The length scale may be equivalent to the correlation 
length scale for fluidity discussed in Refs. [HI, [36[. For 
example, an experimental study of sheared polydisperse 
emulsions found a fluidity length scale comparable to 1-2 
droplet diameters near the glass transition (35| . 

Considering all of our data sets, we do not find a strong 
dependence on either the strain rate 7 or the total strain 
7 for the ranges we consider (7 

mcso 

= 0.006 - 0.02 s~\ 
7moso = 0.3 — 1.6). We do not have a large amount of data 
with which to calculate the correlation functions; unlike 
prior work, we cannot do a time average [73j . If we use 
an exponential function to fit our different data (different 
strains, strain rates, and volume fractions), the mobility 
correlation Sg r yields a length scale ^sr ~ 1-8 — 3.9 /zm 
and the vector correlation Sf yields a length scale w 
3.2 — 7.5 /im. To check this, we also calculate 



(R-S(R)) 

(S(R)) 



(8) 



where the angle brackets indicate averages over R; for 
a perfect exponential, we would have £' = £. Using 
this method, we find more consistent length scales of 
t' Sr « (3.0 - 4.0)a and w (3.5 - 4.3)a. Our data 
do not suggest any dependence of these length scales on 
the control parameters over the range we investigate. Of 
course, as 7 — + 0, we would expect to recover the orig- 
inal unstrained sample behavior 37J. Similar samples 



in this volume fraction range were previously found to 
have length scales with similar values, £$ r « 4a — 8a and 
£,? ~ 6a [73| . However, the time scales for this motion 
are much longer than that for our sheared samples. 



defined as Sf = Ar — (|Af|), in other words, the devia- 
tion of the magnitude of the displacement from the aver- 
age magnitude of all particle displacements. The corre- 
lation functions are computed for the nonaffine displace- 
ments, using At = 45 s to maximize the "signal" (non- 
affine displacements) compared to the "noise" (Brownian 
motion within cages, on a shorter time scale). 

The two correlation functions are shown in Fig. [5] for 
a representative data set. These functions are large at 
short separations R and decay for larger R, suggesting 



C. Defining local plastic deformation 

We wish to look for spatially heterogeneous dynam- 
ics, that is, how the shear-induced motion takes place 
locally and how particles cooperate in their motion. Sev- 
eral prior groups have examined local rearrangements in 
simulations and experiments [HI, [34], HH, [45[, but have 
used differing ways to quantify the motion. We will dis- 
cuss those quantities together, and compare them using 



8 



our data. We have two goals: first, to understand how 
different measures of local rearrangements reveal differ- 
ent aspects of the motion; and second, to see if the spatial 
structure of rearranging groups of particles exhibits any 
particular orientation with respect to the shear direction. 

For all of these definitions of rearranging groups of par- 
ticles, it is useful to define a particle's nearest neighbors. 
Our definition of a particle's nearest neighbors are those 
particles within the cutoff distance tq, set by the first 
minimum of the pair correlation function g(r). 

We start by quantifying the local strain seen by an 
individual particle, which is based on the average mo- 
tion of its neighbors. For a particle i with center at 
fj (t) , the relative positions of its neighbors j are dij (t) = 
r"j(t) — ?i(t). These neighboring particles move, and 
their motions over the next interval At are given by 
Adij (t) = (i + At) - dij (t) , as shown in Fig. [10] 



Adij(t) 
» **». 

\ f 

x ' dij(t+At) 





n(t+At) 



FIG. 10: (Color online) Sketch illustrating the definition of 
the motion of the relative position vector Adij of two neigh- 
boring particles, red and blue (dark and light gray), fi(t) 
represents the position of particle i at time t, dij(t) is the rel- 
ative position between particles i and j at time t, and Ad%j{t) 
is how the vector dij changes over the time interval (t, t + At). 

The strain tensor Ej for this region around particle 
i is then determined by minimizing the mean squared 
difference between the actual relative motions Adij (t) 
and that predicted by Ej, in other words, choosing Ej 
to minimize 



D 



Adij(t) -Eidy(t) 



(9) 



The error, Df ■ , 



quantifies the plastic deformation of 
the neighborhood around particle i, after removing the 
averaged linear response Ejdy(t) [34[. Thus, Df min is 
one way to quantify the local nonaffine rearrangement, 
"local" in the sense of an individual particle i and its 
neighbors. We term D? min the "plastic deformation." 
Note that the sum is computed over the ten nearest par- 
ticles j to particle i, otherwise the value of -D 2 min would 
depend on the number of neighbors. In practice, most of 
these neighboring particles are within 3.0 /im of particle 
i, which is comparable to the first minimum of the pair 



correlation function g(r), which motivates our choice of 
ten neighbors. 

Of course, quite often Ej is different from the overall 
strain over the imaged volume, which in turn is differ- 
ent from the macroscopically applied strain. In practice, 
given that the shear is applied in x direction with the ve- 
locity gradient along z, we only treat the xz components 
of Eqn. 03 that is, dij and Ey can be written as : 



dij — 



,Ei = 



(10) 



To better understand this local strain tensor E,*, we 
follow the method of Ref. [HI, • If the particle-scale 
local strain was identical to the imposed strain, we would 
expect ef z = 7 mcS o and the other matrix elements to be 
zero. We find that these expectations are true on aver- 
age (for example (ef z ) = 7 meS o) but that for individual 
particles their local environment can be quite different. 
For each experiment, the distributions of all four matrix 
elements have similar standard deviations, and examin- 
ing different experiments the standard deviations are be- 
tween 24% - 39% of j meso . 

To quantify the measured particle-scale strain, we de- 
fine the "local strain" 



7- 



%. micro 



(ii) 



(using the definition of the strain tensor which is related 
to E [82]]). That is, this quantity is a local approximation 
to the strain {^f - + t^t)- The local strain 7i, m icro is now 
a second way to quantify the local rearrangement of the 
neighborhood around particle i, in addition to D? min . 
The diagonal elements, ef x and ef z , relate to the dila- 
tion of the local environment. In particular, the local 
environment stretches by a factor of (1 + ef x ) in the x 
direction and likewise (1 + t\ z ) in the z direction. If 
these matrix elements are zero, then the local environ- 
ment remains the same; positive matrix elements corre- 
spond to expansion and negative matrix elements corre- 
spond to contraction. We define the overall dilation as 
Sei — (1 + ef x )(l + ef z ) — 1, which is a third way to 
quantify the local rearrangement around particle i. 

A fourth way to consider local particle motion is the 
previously defined nonaffine displacement, Af 2 . The key 
difference is that -D 2 lin , Tmicroj an d Se all are derived from 
the actual particle motion Ar, whereas Af 2 removes the 
motion caused by 7 meS o (through Eqn. [5]). 

To demonstrate how neighboring particles rearrange 
and result in larger values of these various parameters, 
Fig. [Tj] shows an example using real trajectories. The 
original positions of the particles are shown, along with 
displacement vectors indicating where they move after 
the sample is strained with 7 meS o = 0.58. The overall 
strain is seen in that particles near the top move far- 
ther than those at the bottom; however, the red (dark) 
particle in the middle has an unusual motion, moving 
downward in the —z direction. Figure [TlT b) shows the 
motion of the surrounding particles, as seen in the refer- 
ence frame attached to the red (dark) particle. It is these 
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displacements that are used in the calculation Eqn. [TOl 
Figure [TTT c) shows the predicted final positions of the 
particles (drawn large) based on Ej, as compared to the 
actual final positions (drawn small); the red (dark) par- 
ticle is considered "particle i" . This local region experi- 
ences both shear and a strong dilation in the z direction, 
both captured by Ej. The differences between the pre- 
dicted and actual final positions result in a moderately 
large value of -D 2 lin = 56 /im 2 . In particular, note that 
D 2 lin is defined based on vectors pointing from the red 
(dark) reference particle to the other particles, and be- 
cause the red particle moves downward, the vectors are 
all greatly stretched and this increases Df nin . Finally, 
Fig. fTlT d) shows the positions of the same particles after 
the strain has been applied, where now the box represents 
the mesoscopic strain 7 m0S o = 0.58. The large spheres 
represent the expected positions if the motion was affine, 
and the small spheres show the actual positions. Differ- 
ences between the expected and actual positions result 
in large values of the nonaffme displacement Af. For the 
red (dark) particle, Af = 2.97 /im. Overall, the anoma- 
lous motion of the central particle i is because under the 
large local strain, this particle makes a large jump out of 
its original cage. It is these sorts of unusual motions that 
result in large plastic deformations within the sample. 



D. Collective particle motions 

To investigate the relationships between these quanti- 
ties, a 3 /im-thin y slice of a sample with volume fraction 
4> = 0.51 is shown in several ways in Fig. [T2j In panel 
(a), the x displacement is shown, making the strain ap- 
parent. Panel (b) shows the original Voronoi volumes for 
each particle at t = 0. The Voronoi cell for each parti- 
cle is defined as the volume closer to the center of that 
particle than to any other particle. In subsequent panels, 
the darker colors indicate larger local rearrangements, as 
measured by the non-affine displacement Af 2 [panel (c)], 
plastic deformation £> 2 nin [panel (d)], local strain 7 m i cr o 
[panel (e)], and dilation Se [panel (f)]. It can be seen that 
the darker colored particles cluster together, indicating 
that for each of these measures of local rearrangement, 
the motions are spatially heterogeneous This is a 

real-space picture showing conceptually what is indicated 
by the correlation functions in Fig. [31 that neighboring 
particles have similar motions. These pictures are quali- 
tatively similar to those seen for thermally-induced cage 
rearrangements in supercooled liquids [13, IzS E3| 
and glasses [HI, HH, HH, H3, HI] ■ Furthermore, by compar- 
ing these images, it is apparent that particles often have 
large values of several quantities simultaneously; in par- 
ticular compare panels (c) and (d), and panels (e) and (f). 
While the correspondence is not exact, it suggests that 
all four of these ideas are capturing similar features of re- 
arranging regions. However, it is also clear that there are 
differences between (c) and (d) as compared to (e) and 
(f). In the latter two panels, the region of high activity 




(c) (d) 

FIG. 11: (Color online.) These four sketches show different 
portrayals of a particle with unusual motion, (a) Particle 
motion as seen in the laboratory reference frame. The ar- 
rows indicate displacement vectors, (b) Similar to (a) except 
the motion are in the reference frame where the central red 
(darker) particle is motionless, (c) The large spheres cor- 
respond to the expected positions of the particles, and the 
smaller spheres correspond to the actual positions of the par- 
ticles. The distortion of the box and the expected positions 
of the particles are calculated based on the measured local 
strain tensor E;, where the red (darker) particle is particle 
i. For this local neighborhood, we have 7 m i cr o = 0.73 and a 
dilation primarily in the z direction (e xx — — 0.1, e zz = 0.5, 
and Se = 0.35). (d) The positions of the same particles after 
the strain has been applied, where now the box represents the 
mesoscopic strain 7meso = 0.58. The large spheres represent 
the expected positions if the motion was affine, and the small 
spheres show the actual positions. For all panels, to show 
the displacements in three dimensions better, the radii of the 
large spheres are 0.5 of the real scale. The data correspond 
to a sample with cj> = 0.51, 7meso = 0.58, 7meso = 0.014 s _1 , 
and At = 40 s. 



is spread out over a larger area; more of the particles are 
deforming by these measures at the same time. Nonethe- 
less, in all four cases, the particles around x ~ 2 — 4 
are experiencing the most extreme deformations. 

The Voronoi volume [Fig. [T2T b)] has previously been 
found to be slightly correlated with particle motion in 
unsheared colloidal supercooled liquids [68|, H§|; that is, 
particles with large Voronoi volumes have more space to 
move and thus are likely to have larger displacements 
than average. Here it appears that there is no correla- 
tion between the Voronoi volume and the particle motion, 
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(a) Ax (b) Voronoi 



z 




FIG. 12: A 3 /im thin cut through the sample, showing the 
xz shear plane; the axes are as labeled in (a). The particles 
are drawn in their location at t = 0, the start of this strain 
half-cycle, (a) Darker particles have larger values of Aa:, thus 
indicating the applied strain. The top particles (large z) are 
moving left (dark colors) and the bottom particles are moving 
right (light colors), (b) Darker particles have larger Voronoi 
volumes at t = 0. (c) Darker particles have larger values of the 
non-affine motion Af 2 . (d) Darker particles have larger values 
of the plastic deformation D^ in . (e) Darker particles have 
larger values of the local strain 7 m icro- (f) Darker particles 
have larger values of the dilation Se. The sample is the same 
as the data shown in Fig. O see that figure caption for details. 



suggesting that for these strained samples the local vol- 
ume is not a crucial influence on the motion [38[. This 
is probably because in the end, all cages must rearrange 
for the strain to occur. 

To demonstrate the relationships between the mea- 
sures of plastic deformation more quantitatively, we com- 
pute 2D histograms comparing pairs of the variables, 
shown in Fig. [T5] The darker color indicates larger joint 
probability, and the dotted line represents the mean value 
of the quantity on the vertical axis corresponding to the 
quantity on the horizontal axis. Figure [T3T a) shows that 
on average, particles with a large plastic deformation 



D^ in are also much likelier to have a large nonaffine dis- 
placement Af 2 . This is suggested by the specific example 
shown in Fig. [TTTc.d). Similarly, Fig. [TBT b) shows that a 
particle's microscopic strain 7 m i C ro is well correlated with 
the dilation Se. For these data, the mesoscopic strain is 
7mcso = 0.43; particles with 7 micro < 7 moso are more of- 
ten in local environments that contract (Se < 0), and 
vice-versa. As a contrast, Fig. fT3Tc) shows a somewhat 
weaker correlation between -D„ lin and 7 m i cr o- The Pear- 
son correlation coefficients between these quantities are 
C(D 2 min , Af 2 ) = 0.43 ± 0.11, C( 7micro , Se) = 0.42 ± 0.09, 
and C(7 m i cro , D^ in ) = 0.17 ±0.05. The uncertainties are 
from the standard deviations of the correlation coeffi- 
cients from the nine different experiments we conducted. 

Overall, Fig. [T31 suggests that -D 2 ^ and Af 2 both cap- 
ture the idea of plastic deformation [661 ] . The correspon- 
dence between these two variables is nontrivial, given 
that Af 2 is based on trajectories with the overall strain 
removed (the strain computed from all observed parti- 
cles), whereas -D 2 lin only accounts for the very localized 
strain of the neighboring particles. In contrast, 7 m icro 
and Se are well-suited to examine particles moving in 
atypical ways; typical particles have 7 m icro ~ 7mcso and 
Se = 0. These two separate ideas (plastic deformation 
and atypicality) are only weakly correlated. Other than 
the three specific comparisons shown in Fig. [T3] all other 
comparisons are even less correlated (|C < 0.1). We also 
examined the quantities |7 m i cr o — 7meso| and \Se\ as ways 
to measure the deviations from typical behavior; these 
quantities are also only weakly correlated to the other 
measures of deformation. 

We now return to the question of the isotropy of the 
motion. Figure [H] indicates that the distribution of all 
particle motions is isotropic, but it is possible that the 
spatially heterogeneous groups of highly mobile particles 
shown in Fig. [12] are themselves oriented along a pre- 
ferred direction. To investigate the 3D structures of these 
relaxation regions, we quantify the sizes of these active 
regions in the x, y and z directions. To start, we define 
connected neighboring particles as those with separations 
less than ro, the distance where pair correlation function 
g(r) reaches its first minimum. (Note that this is slightly 
different from the neighbor definition used for Eqn.[9l see 
the discussion following that equation.) For a given quan- 
tity, we consider active particles as those in the top 20% 
of that quantity, similar to prior work [3?], 5H, [ZZl [85| . We 
then define the active region as a cluster of connected ac- 
tive particles. For example, Fig. [TJ] shows a cluster of 
particles with large non-affine displacements [panel (a)] 
and a cluster with large plastic deformations [panel (b)] . 
Each cluster is drawn from the same data set, and the 
particles drawn in red (darker) are common to both clus- 
ters. (Note that clusters drawn based on 7 m i cr o and Se 
are smaller. In Fig. H"27 e,f), more regions have large val- 
ues of these parameters, but the top 20% most active are 
not clustered to the extent they are in Fig. [T2Tc.d).) 

We wish to understand if the shapes of such clusters 
show a bias along any particular direction. It is impor- 
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(a) Ar 2 cluster 




0.5 1 

Ymicro 

FIG. 13: (Color online) 2D histograms of (a) D^ in and Af 2 , 
(b) 7micro vs. dilation 8e and (c) 7 m icro and D„iin- The 
sample has the same parameters as Fig. [5] see that caption 
for details. In these histograms, darker colors stand for the 
larger probability. The dotted lines are the average values of 
the y axis corresponding to each value on the x axis. The 
correlation coefficients between each pair of variables are (a) 
C = 0.43±0.11, (b) C = 0.42±0.09, and (c) C = 0.17±0.05; 
see text for details. 



(b) Dmin cluster 

FIG. 14: (Color online.) A cluster of neighboring particles 
with large values of (a) non-affine motion Af 2 and (b) plastic 
deformation D^ n . To avoid biasing the cluster shape, the 
particles are drawn from within a 15x15x15 pm 3 cube (solid 
lines). The dashed lines indicate the spatial extent of the 
two cluster shapes in each direction. The cluster in (a) is 
14 x 12 x 11 /im 3 and the cluster in (b) is 8 x 14 x 10 /mi 3 . 
The red particles are common between the two clusters, and 
the black bonds indicate neighboring particles; particle sizes 
are drawn 0.8 times the actual scale, to make the connections 
more visible. In each case, the particles shown are the top 
20% for the parameter chosen. The data corresponds to the 
same parameters as in Fig. [5] see that caption for details. 



tant that the experimental observation volume not bias 
the result, so from within the observed 3D volume we 
consider only particles that start within an isotropic cube 
of dimensions 15x15x15 /xm 3 . The size of this cube is 
chosen to be the largest cube for which all the particles 
are within the optical field of view for the full half-cycle 
of shear for all experiments. (See the related discussion 
at the end of Sec. [TT| ) Within this isotropic volume, we 
consider the largest cluster from each experiment and for 
each considered variable. We then define the extent of 
that cluster in each direction from the positions of each 
particle within the cluster as ^extent = max(cci) — mm(xi) 
and similarly for y and z. 

Anisotropic cluster shapes would be manifested by sys- 



tematic differences in the relative magnitudes of ^extent; 
J/extent, and z cx tcnt- We compare these in Fig. [15] for clus- 
ters of particles with large non-affine motion Af 2 [panel 
(a)] and large plastic deformation -D 2 ^ [panel (b)]. The 
comparison is made by using the ratios y C xtcnt /^extent 
and Extent /^extent , thus normalizing the extent in the 
y and z directions by that of the shear velocity direc- 
tion x. Thus, if a cluster has the same extent in x and 
y, J/extent /^extent should be equal to 1, along the ver- 
tical dashed line. Similarly, for the same extent in x 
and z, points should be along the horizontal dashed line 
with Zextent /^extent = 1- If the extent is the same for y 
and z, the points should be along the diagonal line with 
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2/oxtcnt /xoxtout = Extent /^extent- For an isotropic cluster 
with same size in all dimensions, the point should be in 
the center (1,1). 

As shown in Fig. 1151 for all of our data, we find no sys- 
tematic anisotropy; the cluster extent ratios are mostly 
clustered around the isotropic point (1,1). Due to ran- 
dom fluctuations, no cluster is perfectly isotropic, yet the 
points seem fairly evenly distributed around the three 
dashed lines. Thus, while the shear-induced rearrange- 
ments take place in localized regions (Fig. [T2")) . the data 
indicate that these regions on average have no directional 
bias. This seems true for both the nonaffine displace- 
ments Af 2 in Fig.fToTa) and the plastic deformation -DjL, 
in Fig.H^b). 
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FIG. 15: Comparison between the cluster extent in x, y and 
z, based on different variables, (a) Particles with the largest 
non-affine motion Af. (b) Particles with the largest plas- 
tic deformation D^ in . The symbols represent different vol- 
ume fractions: <j> = 0.51 (circles), 0.56 (triangles) and 0.57 
(squares). Colors (or grays from light to dark) represent dif- 
ferent accumulative strains 7meso = 0.3 — 0.9. The clusters 
are comprised of the top 20% of the particles with the given 
characteristic. 



IV. DISCUSSION 

We examined the microscopic plastic deformations oc- 
curring in several sheared dense colloidal suspensions. 
Our first main observation is that on average, individ- 
ual particles have no bias in their direction of motion, 
other than that trivially imposed by the strain. When 



this imposed motion is removed from the particle trajec- 
tories, the remaining shear-induced motion is isotropic: 
particles are equally likely to move in any direction. Our 
second main observation is on the shape of groups of par- 
ticles undergoing plastic rearrangements. There are sev- 
eral ways to determine which particles are rearranging, 
and we have shown that all of these are useful for high- 
lighting local regions of deformation. Furthermore, the 
shapes of these regions are also isotropic. However, we 
cannot rule out that with more data and subtler analysis, 
we might find anisotropies in particle motion 90]. 

In our results, we find little dependence on the overall 
volume fraction 0, total strain, or strain rate. For the 
volume fraction, all of our samples are dense liquids with 
4> < 4>g- At significantly lower volume fractions, presum- 
ably particles would not be caged and the shear-induced 
rearrangements might be quite different [70l ]. At higher 
volume fractions 4> > fid prior work has seen similar re- 
sults [To|, [39( although not examined the shapes of the 
rearranging regions in detail. It is possible that results 
in glassy samples might be different, given that near <Pg 
slight changes in volume fraction have large consequences 
[42j |. but we have not seen clear evidence of that in our 
data. For the total strain, we have not examined a wide 
range of parameters. In all cases, we are studying suffi- 
ciently large enough strains to induce irreversible, plas- 
tic rearrangements. For the strain rates, all of our strain 
rates are fast enough such that the modified Peclet num- 
ber Pe* is at least 7, so that thermally induced diffusive 
motion is less relevant. It is likely that at slower strain 
rates (lower Peclet numbers) , different behavior would be 
seen [37| . 

Previous work [54], |9j H^, [9^| found that oscillatory 
shear can induce crystallization of concentrated colloidal 
suspensions. The 'induction time' of this crystallization 
is strain dependent: a larger strain amplitude results in 
shorter induction time. In our experiments, we studied 
only a limited number of oscillations, and our strain am- 
plitude ~1. We did not observe crystallization in any of 
our experiments. It is likely that were we to continue our 
observations for much longer times, we could see the on- 
set of shear-induced crystallization, and so we note that 
our experiments are probably studying a non-equilibrium 
state. Additionally, Fig. 2] shows that our strain rate 
takes a while to stabilize after flow reversal, which again 
suggests that our results are not in steady-state. Thus, it 
is possible that our primary observation, that the shear- 
induced particle rearrangements are isotropic in charac- 
ter, is limited only to the transient regime we observe. 
It is still intriguing that in this regime, particle mo- 
tion is so isotropic. For example, Fig. [4] shows that the 
sample takes a while to requilibrate after shear reversal, 
yet there is no obvious signature of this in the particle 
motion or the configurations of the particles. Likewise, 
presumably the long-term crystallization will be caused 
by anisotropic motion (and result in further anisotropic 
motion), but no signs of this are present in the early- 
time amorphous samples we study. It would be inter- 
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esting to conduct longer-term experiments to relate the 
particle rearrangements to those resulting in crystalliza- 
tion. Alternatively, it would be also interesting to use a 
cone-and-plate geometry shear cell capable of indefinitely 
large strains [62| , to reach the steady-state constant shear 
regime. 
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